Two loop lattice expansion of the Schrodinger 
functional coupling in improved QCD 



A. Bode P. Weisz ^ and U. Wolff " 

"^SCRI, Florida State University, 
Tallahassee, Fl, 32306-4130, USA 

''Max-Planck-Institut fiir Physik, Fohringerring 6, 
D-80805 Miinchen, Germany 

'^Institut fiir Physik, Humboldt Universitat, Invalidenstr. 110, 
D-10099 Berlin, Germany 

February 1, 2008 



Abstract 

The contributions of the improved fermion action of Sheikholeslami 
and Wohlert to the two loop coefficient of the expansion of the Schrodinger 
functional coupling in terms of the lattice coupling are calculated for the 
gauge group SU(3). These coefficients are required for the second order 
relation of lattice data to the MS-coupling. By taking into account all im- 
provement coefficients we are able to improve the Schrodinger functional 
to two loop order. This article is based on the revised version of 0|. 



1 Introduction 

In the framework of the Schrodinger functional the renormalized coupling asF 
can be traced from low to high energy numerically on the lattice with finite size 
techniques. At high energy a conversion to with perturbation theory is 
possible. The relation can be calculated by expanding both and asp in 9o, 
the lattice bare coupling. 

This program has been completed for the pure gauge sector (A^f = 0) of QCD 
where the two loop relation between and asF was required in order to 
avoid a significant source of error in the conversion to ctj^. Here we present 
the perturbative two loop relation for arbitrary N{ for the mass-independent 
scheme described in ||^. 
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The outline of the calculation is analogous to the pure gauge case We 
have to calculate the perturbative coefficients depending on the box size / = 
L/a, where asF is defined. From this we are able to extract the continuum 
results and the 0(a) improvement terms. 



2 Definition of a 



SF 



The coupling is defined via the effective action F in a finite box with extension 
L 

exp(-F) = j D[U]D{ii\D[^]eM-S[U,'^M (1) 

where S consists of the Wilson gauge action and Shcikholcslami-Wohlcrt fcrmion 
action. In the spacelike directions periodic boundary conditions with an addi- 
tional phase for the fcrmions are used. The timelike boundaries are given by 
a constant and diagonal colour field depending on the parameter 77. Certain 
1 ±7o projected fermion field components are set to zero at t = 0, L. To achieve 
0(a) improvement we have to add coupling dependent weights to some links 
and plaquettes on the boundaries. Details of the definition of the boundaries 
can be found in Q . 

The coupling asF(g) = gl^{L) / {A-k) , at g = 



F'(i,r7) 



(2) 

r)=0 



is normalised via the classical action minimal Fq. The structure of the expansion 
is given by 

9h{L) = gl + bio(/) +pii(/)iVf]5o' (3) 
+ b2o(/) +P2i{I)Ni+p22{I)Nf]gl + 0{gl) 

where the dependence on the one and two loop improvement coefficients is 
hidden in pij. The pij are sums of Feynman diagrams. The coefficients presented 
here arep2i(^) and P22{I)- The one loop coefficient pii{I) was computed in |^. 



3 Perturbative expansion 

Since the classical action minimum is not affected by the presence of the fermions 
we can use the gluon and ghost propagators of the pure gauge case. The fermion 
propagators must be calculated numerically. A recurrence scheme similar to 
the other propagators is developed and will be described elsewhere The 9 
diagrams involving fermions are summed in position space. 

Several discrete symmetries can be used to reduce the number of terms. 
The computational complexity and resulting memory requirements increase by 
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including the fermions due to the Dirac trace. Hence we coded two versions, 
one where we store all components of the Dirac propagator and traces are easy 
to evaluate and one with a reduced set of spinor components (due to the cubic 
symmetry in space) but slightly more complicated Dirac traces. 

Careful tests were applied: the symmetries of propagators, diagrams and re- 
sults, gauge parameter independence of p2j{I), comparison of results for small 
L with an independently written program using momentum space sums, verifi- 
cation of the universal /3-function coefficients, consistency relations of Symanzik 
improvement. 

In order to get the continuum relation we have to extract various (sub)leading 
coefficients of the asymptotic series in a/L describing p2i- Apart from the well 
known blocking procedure an alternative method has been developed and will 
be described in Since further development is in progress and further data 
are being generated, all results are preliminary. 

This gives up to 0{\n^{I)/P): 

Pw{I) = 26oo ln(/) + 0.36828215(13) , (4) 

pii(/) = 26oi ln(/) - 0.034664940(4) , (5) 

P2o{I) = 2hw HI) + 0.048085(63) + pio{lf , (6) 

P2i{I) = 26ii ln(/) - 0.004738(56) + 2pio(/)pii(/) , (7) 

P22(/) - 0.00 0211(29) -fpii(/)2 . (8) 

All terms oc ln*(/)// are zero due to improvement. The fcrmionic contributions 

(2) 

to the boundary improvement coefficients c\ are also determined. In addition 
the expansion coefRcients for the standard Wilson action Cgw = 0, are calculated. 
This is numerically cheap because the clover vertices dominate the calculation. 



4 Evolution of agp 

The non-perturbative evolution of the coupling agp for A^f = 0| and the 
perturbative three loop evolution for A'f = 0,2 is shown in figure 1. These 
curves are integrated with the smallest asF value as initial value (the absolute 
scale applies only for A^f = 0). From perturbation theory we can expect a visible 
shift in the evolution of the coupling going from A'f = to A't = 2. 



3 



0.3 
0.25 
0.2 

asF(q) 

0.15 
0.1 
0.05 

1 10 100 

q [GeV] 

Figure 1: Evolution of the coupling asF- 

5 Connections to other couplings 

Up to now Oj^jg in ckq is known to two loop order only for the unimproved action, 





nonpert. n,=0 i — i — i 




fit to nonpert. 




Three Loop, nf=0 




Three Loop, nf=2 







IBI, |9||. Connecting the series for N{ = 2 gives 

"ms(^9) = "SF + Cl(s)asF('?)^+C2(s)QsF(<?)^ , (9) 

ci(s) = -87r6oln(s) + 1.3353 , (10) 
C2(s) = ci(s)2 ~ 32TT^bi ln(s) + 3.13(3) . (11) 

With known we are able to quote 

bf{AT:f = 0.482(7) - 0.275(8)iVf + 0m6{2)Nl - 0.00175(8)iVf^ (12) 

which results in = 0.063(24)/(47r)3 for iVf = 2. The perturbative relation 
(^ between asF and involves the scale factor s. Fixing it by demanding 
either ci(s) = or C2(s) to be minimum gives comparable results: 



a^(2.3820g) = asF(g) + 0.50(3)asF(g)' , 

a^(2.9243g) = asF(g) - 0.3156asF(9)' + 0.40(3)asF(g)' 



(13) 

(14) 



where the quantities quoted without error are in all digits significant. The large 
two loop coefficients might require smaller couplings for a conversion to obtain 
the same error as for the iVf = case. 



4 



The two loop relation between the lattice couplings for the improved theory 
and the Wilson fermion action and the corresponding three loop /3-function 
coefficients are (for N = 3) also available. 



6 Finite a effects 



In numerical simulation the finite a step scaling function 

9sf{sL) = T,{s,glp{L),a/L) 
is accessible, and the continuum limit 

(^{s,9sf{L)) = lim'E{s,gy{L),a/L) 

a— »0 



(15) 



(16) 



has to be taken. The perturbative finite a effects are defined from the series pij 
via 



(17) 



The ^'s hopefully give a hint about the magnitude of scaling violation in nu- 
merical simulations. Scaling for a in perturbation theory looks comparable for 
N{ = 2 and A^f = 0. This is shown in figure 2 for the two loop contribution 
62{2,a/L). 



LVa 



inf.16 119 8 7 




-2 
-4 
-6 

o 

^ -8 
_i 

"a 

cvT -10 
to 

-12 
-14 
-16 
-18 



-1 1 1 — I r- 



Nf=2, full improvment + 

N(=0, full improvment x 
_i ■ ■ I ■ ■ I I ■ ■ ■ I ■ b_ 



0.01 



0.02 0.03 0.04 0.05 0.06 0.07 

,2 



Figure 2: Two loop contribution to the scaling 
violations of the step scaling function a. 
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7 Conclusion 



The two loop expansion of the Schrodinger functional with the Sheikholeslami- 
Wohlert fermion action is numerically demanding in comparison to the pure 
gauge case. For agp the three loop /3-function is now known and conversion to 
a-j^ is possible at two loops. For the dynamical simulations of the Schrodinger 
functional the perturbative results have small 0(a) effects and do not indicate 
particular complications for the continuum extrapolation. 
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